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Abstract 



I> ' Oki Island is located between Japan and Korea along the Tokai-To-Kamioka 

(T2K) baseline. The distance from J-PARC to Oki Island is about 653km, which is 
twice that of the T2K experiment (L = 295km). When the off-axis angle of the neu- 
trino beam from J-PARC is 3.0° (2.0°) at Super-Kamiokande (SK), the off-axis beam 



On 

(OAB) with 1.4° (0.6°) reaches at Oki Island. We examine physics case of placing 
a far detector in Oki Island during the T2K experimental period. We estimate the 
matter density profile along the Tokai-to-Oki baseline by using recent seismological 
measurements. For a detector of 100 kton fiducial volume and 2.5 x 10 21 POT (pro- 



tons on target) exposure for both and beams, we find that the mass hierarchy 
pattern can be distinguished at 3 a level if sin 2 2# RCT = 4|£/ e 3| 2 (l — |t/ e 3| 2 ) <; 0.09, by 
observing the electron-like CCQE (Charged-Current Quasi Elastic) events. The CP 
phase in the Maki-Nakagawa-Sakata lepton flavor mixing matrix, <5 MNS , can be con- 
strained with ±20°. As a reference, we repeat the same analysis by placing the same 
detector in Korea at L = 1000 km and OAB=0.5° (T2KK) and also by placing it at 
the SK site (T2Ki22)- The Tokai-to-Kamioka-OKI (T2KO) sensitivity to the mass 



hierarchy is about 1/3 (in ^Xmin) °^ T2KK, while the sensitivity to the phase 5 



MNS 



is similar between T2KO and T2KK. The T2K122 option has almost no sensitivity 
to the mass hierarchy, and cannot measure the CP phase except when 5 MNS ~ —90° 
(90°) for the normal (inverted) hierarchy. 
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1 Introduction 



Since the neutrino oscillation was first observed at Super-Kamiokande (SK) p], many 
experiments have measured the physics parameters of the neutrino oscillation [2]-|12j. 

Some of these experiments observe the survival probability of and which are gener- 
ated in the atmosphere by the cosmic ray [2J. Accelerator based long baseline experiments 
[3llll|5] also measure the survival probability. From the combined results of these 
experiments [2J-[5], the mass-squared difference and the mixing angle are obtained as 

sin 2 2fl ATM > 0.90 (90%C.L.) , (la) 
\5m 2 ATM \ = (2.35i°^) x 10- 3 eV 2 . (lb) 

The sign of the mass-squared difference, eq. (llbj) . cannot be determined from these exper- 
iments. The SK collaboration also reported that the atmospheric neutrinos oscillate into 
active neutrinos [6]. 

The combined results of the solar neutrino observations [3 [8], which measure the survival 
probability of v e from the sun, and the KamLAND experiment [9], which measure the v e 
flux from the reactors at distances of a few 100 km, find 

sin 2 2# SOL = 0.852l°;°26> (2a) 
Kol = (7.501^) xlO" 5 eV 2 , (2b) 

where the sign of mass-squared difference has been determined by the matter effect inside 
the sun [13]. The SNO experiment determined that v e from the sun changes into active 
neutrinos [8]. 

For the reactor experiments, which measure the survival probability of v e from the 
reactor at distances of L ~ 0(1) km [TUl [TT] . the CHOOZ experiment observed no signal, 
giving 

sin 2 2# RCT < 0.17, (3a) 
for 1 5ml CT | = 2.35 x 10" 3 eV 2 , (3b) 

at the 90% confidence level [TU]. More recently, the T2K collaboration reported several 
Vp, —> v e appearance candidate events giving 

0.03 (0.04) < sin 2 26 RCT < 0.28 (0.34) (4) 

for sin 2 # 23 = 0.5 and S MNS = 0° with the normal (inverted) hierarchy at the 90% confidence 
level [T4]. The MINOS collaboration reported 

sin 2 26 RCT = 0.041^(0.079^1) (5) 

also for sin 2 ^23 = 0.5 and S MNS = 0° with the normal (inverted) hierarchy [15]. The Double 
CHOOZ collaboration, which is one of the new reactor experiments, found hints of reactor 
electron anti-neutrino disappearance consistent with neutrino oscillation and reported 

sin 2 26 RCT = 0.086 ± 0.041 (stat.) ± 0.030 (syst.) (6) 



from analyzing both the rate of the prompt positrons and their energy spectrum [T6] . 
Recently another new reactor experiment, the DayaBay experiment [17], announced that 
they have measured the neutrino mixing angle as 



sin 2 29 RCT = 0.092 ± 0.016 (stat.) ± 0.005 (syst.) 



(7) 



which is more than 5a away from zero. The RENO collaboration, which also measure the 
reactor v e survival probability, shows the evidence of the non-zero mixing angle; 

sin 2 20 RCT = 0.113 ± 0.013 (stat.) ± 0.019 (syst.) , (8) 

from a rate-only analysis, which is 4.9cr away from zero. 

Since the MiniBooNE experiment [12] did not confirm the LSND observation of rapid 
— > v e oscillation [18], there is no clear indication of experimental data which suggests 

and 



more than three neutrinos. Therefore the — > v e appearance analysis of T2K 
MINOS pj)] presented above assume the 3 neutrino model, with the 3x3 flavor mixing, 
the MNS (Maki-Nakagawa-Sakata) matrix [19] 



(9) 



relating the weak interaction eigenstates (z/ e , v^, v T ) and the mass eigenstate Pi with the 
mass rrii (i = 1,2,3). The mass-squared differences that dictate the neutrino oscillation 
phase are then identified as 
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(10a) 
(10b) 



where only the magnitude of the larger mass-squared difference, m\ — m\ is determined 
in eq. ( 11b|) . The m\ > m\ case is called normal, while m\ < ml case is called inverted 
neutrino mass hierarchy. With a good approximation [20], the three independent mixing 
angles of the MNS matrix can be related to the oscillation amplitudes in eqs. ( ITa| . ff2al) . 
and sin 2 29 RCT in eqs. (|3al) . (Hj) to ((71), with the MNS matrix elements U a i : 



sin 9 A tm 
sin 6» RCT 
sin 29 SOh 



U, 3 

\U e3 \ 

2U el U e2 



sin 9 23 cos 6>i 3 , 
sin 6»i3 , 

sin 2# 12 cos 2 ^13 



(11a) 
(lib) 
(11c) 



In the last equations, the defining region of the three mixing angles 9^ = 9ji can be chosen 
as < #i2, 6*i3, 6*23 < vr/2 [21], which is consistent with the convention of non-negative U e \ 
and U e 2 [22] ■ The argument of U e3 gives the CP phase of the MNS matrix, 



5mns = -arg[/ e3 



(12) 



Even after 9 RCT , the smallest of the three mixing angles, is measured by the accelerator 
based long baseline experiments [TH [15] and the reactor experiments [HI [T7J [23], three 
unknowns remain in the 3 neutrino model, which are the sign of the larger mass-squared 
differences, normal (ml — ml > 0) or inverted (m 2 — < 0), the leptonic CP phase (<5 MNS ), 
and the sign of sin 2 9 ATM — 0.5 if its magnitude differs from zero significantly. The main 
purpose of the next generation neutrino oscillation experiments is to determine these three 
unknown parameters of the three neutrino model. 

The accelerator based long baseline neutrino oscillation experiment with two-detectors 
for one-beam, such as Tokai-to-Kamioka-and-Korea (T2KK) experiment [21]- [30], is one of 
the promising experiments for measuring all the three unknowns. When one measures the 
neutrino energy (E v ) and the magnitude of the v e appearance probability at significantly 
different baseline lengths, the degeneracy between the sign of the larger mass-squared 
difference and the CP phase can be resolved [26]- [30], since they affect the magnitude of 
the v e appearance probability and the neutrino energy at the first oscillation maximum 
differently at the two baseline lengths. Once the mass hierarchy and the CP phase are 
determined, the sign of sin 9 ATM — 0.5 can be determined |28] since the — y v e oscillation 
probability is proportional to sin 2 9 ATM sin 2 9 RCT . 

In the previous works [26]- [30], physics impacts of the T2KK experiment have been 
studied systematically, and the following observations have been made: If a 100 kton 
water Cerenkov detector is placed in Korea at L = 1000 km observing the 0.5° off-axis 
beam (OAB) during the T2K exposure time of 5 x 10 21 POT (protons on target), the 
mass hierarchy can be resolved at 3cr level for sin 2 29 RCT ;> 0.05 (0.06) when the hierarchy 
is normal (inverted), and the CP phase can be constrained uniquely, by measuring the 
CCQE (Charged- Current Quasi Elastic) events [26J [27]. The sign of sin 2 9 ATM — 0.5 can 
also be determined for | sin 2 9 ATM — 0.5| = 0.1 with 3a level, if sin 2 29 RCT > 0.12 for the 
normal hierarchy [28] . When we take into account the smearing of reconstructed neutrino 
energy due to finite detector resolution and the Fermi motion of target nucleons, resonance 
production, and the neutral current (NC) n° production background to the v e appearance 
signal, it is found that the mass hierarchy pattern can still be determined at 3a level for 
sin 2 2# RCT ;>0.08 (0.09) when the hierarchy is normal (inverted), but the CP phase can no 
longer be established at 3a level with 5 x 10 21 POT exposure [29]. In Ref.[30], matter 
distribution profile along the T2K and Tokai-to-Korea baselines have been studied, and 
merits of splitting the total exposure time into half and half beam have been reported. 
Studies on the impacts of the neutrino energy smearing and the NC ir° background for the 
v e oscillation measurements are in progress [3Tj . 

Oki Island^ is placed between Japan and Korea along the T2K baseline and inside 
Japanese territory. The distance from J- PARC to Oki Island is about L = 653 km which is 
about two times longer than that of the T2K experiment [L = 295 km). In Fig. [lj we show 
the surface map of the T2K, T2KO (Tokai-to-Kamioka-and-Oki), and T2KK experiments, 

*This island is sometimes called "Oki-no-Shima" , because "Shima" means an "Island", "no" is "in", 
and "Oki" means the "Offing" ; an island in the offing, or an offshore island. 




Figure 1: The surface map of the T2K, T2K0, and T2KK experiment. The yellow blobs 
show the center of the neutrino beam for the T2K experiment at the sea level, where the 
number in the white box is the off-axis angle at SK. 

in which the red lines denote the baselines for each experiment and the yellow blobs show 
the center of the neutrino beam from J-PARC at the sea level, when the off-axis angle 
at SK is 2°, 2.5°, and 3°. The authors of Ref.[3"2"] consider the physics performance of 
a 100 kton Liquid Argon Time Projection Chamber placed at Oki Island. In Ref. [33], 
the authors studied the Oki Island site from the geological, geographic and infrastructure 
points of view for the possibility to construct a large detector. They conclude that Oki 
Island is one of the good candidate sites for a large detector. 

In this paper, we study the physics potential of a 100 kton water Cerenkov detector 
placed in Oki Island instead of Korea during the T2K experimental period, by using exactly 
the same setting assumed for the T2KK experiment in Ref. [30J, except for the location of 
the detector. This allows us to compare the physics capability of the two proposals on the 
same footing. Since we do assume in this analysis that the — > v e CCQE events can 
be distinguished from the background, our results should equally be applied for a more 
advanced Liquid Argon detector [32] . 

This article is organized as follows. In the next section, we briefly review neutrino 
oscillation in the matter and give useful approximation formula for — > z/ M , v e and — > 
z/ M , v e oscillation probabilities. In section [31 we show the matter profile between J-PARC 
and Oki Island by using recent seismological measurements and give relations between the 
off-axis angle observable at Oki and that at SK. The expected numbers of — > z/^z/g 
[y^ — > z/ M , v e ) CCQE events at the far (Oki) and the near (SK) detectors are shown for 
typical parameters in section HJ In section [5j we present the x 2 function with which 
we estimate the statistical sensitivity of the T2KO and T2KK experiments on the mass 
hierarchy and the CP phase. In section [6j we show the results on the mass hierarchy 
determination of the T2KO experiment and compare it with T2KK, and also with T2K122 
where the 100 kton detector is placed in the Kamioka site making the total fiducial volume 
122 kton at L = 295 km. In section [3, we show the CP phase sensitivity of T2KO, T2KK, 



and T2K122 experiments. We summarize our findings in the last section. 



2 Neutrino oscillation in the matter 



In long baseline neutrino oscillation experiments, neutrinos from the accelerator interact 
coherently with electrons and nucleons by charged and neutral current interactions. These 
coherent interactions make an additional potential in the effective Hamiltonian [13J. Be- 
cause the potential from the neutral current interactions are flavor-blind, it does not affect 
the neutrino oscillation probability in the three neutrino model. On the other hand, the 
charged- current interactions contribute only to the effective potential of u e and u e . The 
neutrino oscillation probabilities hence depend on the matter density profile along the base- 
line. Detailed studies on the matter distribution dependence of the oscillation probabilities 
for the T2K and T2KK experiments are found in Ref. [30]. In this exploratory analysis, 
we compare the physics capability of the T2KO, T2KK, and T2K experiments by using 
the average matter density along the baselines which is found to approximate the matter 
effects rather accurately [30] . 

The Hamiltonian of a neutrino propagating in the matter is then expressed as 
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on the flavor space (u e , u^, u T ) T , with 5m? 
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(13a) 



(13b) 
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(14) 



gives the matter effect with the electron number density n e , which is approximated by the 
average matter density p (g/cm 3 ) along the baseline. By using the solution of eq. (I13bj) . the 
oscillation probability that an initial flavor eigenstate u a is observed as a flavor eigenstate 
up after traveling a distance L along the baseline is 
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where 



A- = A ' 



2E„ 



(16) 



All our numerical results are based on eq. ( !T5|) . However, we find the following analytic 
approximations [26J [30] useful for understanding the reason why and how the one-beam two- 
detectors long baseline experiments such as T2KK |26j- [30] and T2KO, with a far detector 
in Oki Island, can determine the neutrino mass hierarchy and the CP phase simultaneously. 
They are obtained by expanding the oscillation probabilities in terms of the three small 
parameters; the matter effects, which is small at energies below a few GeV around the earth 
crust for p ~ 3 g/cm 3 and L <^ 1000 km, the oscillation phase A 12 = 5m\ 2 L/2E v , which is 
also small (~ 7r/30) near the first oscillation maximum, |A 13 | ~ it, and the mixing factor 
\U e 3\ 2 = sin 2 # RCT , eq. ( 11 lb}) , which is about 1/40 from the recent reactor measurements 
eqs. © and (ED. 

The survival probability can then be approximated as 

iV", = 1 - sin2 ^atm (1 + A») sin 2 f ^±^ \ , (17) 
around the first maximum of the main oscillation phase 

A 13 = 5m 2 = — -L , (18) 

3 13 2E„ 2E tJ ' v ; 



with 



A» = -^-(l-tan 2 ^ ATM )sin 2 ^ CT , (19a) 

B» = ^(l-tan 2 # ATM )sin 2 fl RCT (19b) 
-A 12 (cos 2 6 SOL + tan 2 

atm sin 2 ^sol sin 2 R ct — tan ^ ATM sin 29 SOL sin 9 KCT cos 5 mns ) 

We find that the above formula reproduce the survival probability with 1% accuracy 
throughout the parameter range explored in this analysis, which covers all the three neu- 
trino model parameters in the 3cr allowed range of the present measurements eqs. (JIJ-flH]), 
for the neutrino energies 400 MeV < E v < 4 GeV (for 9^), and for the baseline lengths 
295 km < L < 1200 km, except where the probability is very small, (P v _>.„ <, 10~ 5 ) 
, [30] . The matter effects are proportional to the term 

^-■OJBaf')!-^-) (20) 



2E V V.3g/ cm ) V 1000km / 

and it appears with the opposite sign both in the amplitude correction term eq. (119a}) 
where A 13 = |A 13 | for the normal while A 13 = — |A 13 | for the inverted hierarchy, and also 
in the phase shift in eq. ( ITTj) . where |Ai 3 + B^\ = n + B^ for the normal, tc — B 11 for the 
inverted hierarchy around the oscillation maximum |Ai 3 | ~ it. Nevertheless, the effect is 
too small to be observed in the near future because the term 1 — tan 2 9 ATM is bounded as 

0.48 > 1 - tan 2 9 ATM > -0.92 (21) 



at 90% C.L. from eq. ( ITaj) . and because both A 11 and B^ are proportional to the product of 
the two small parameters, the matter effect term of a§Lj2E v and sin 2 9 RCT . The coefficient 



of cos5 MNS in is also small, only 1% of the main oscillation phase |Ai3|. In other words, 
the survival probability is very insensitive to the unknown parameters of the three 
neutrino model, and hence can be used to measure sin 2 $atm and | 5?ti 2 3 | accurately 

Under the same conditions that give eq. (117j) for the survival probability, the u e 
appearance probability can be approximated as [30] 
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where we retain both linear and quadratic terms of A 12 and do- The analytic expressions 
for the correction terms A e , B e and C e are found in Ref.[30j. For our semi- quantitative 
discussion below, the following numerical estimates [30] for sin 2 29 ATM = 1 and sin 2 2# s 
0.852 suffice: 
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The term C e is relevant only when the — >■ v e oscillation probability is very small [30J. 

The first term in A e in eq. (I23aj) is sensitive not only to the matter effect but also 
to the mass hierarchy pattern, since Ai 3 ~ n for the normal while Ai 3 ~ — tt for the 
inverted hierarchy around the oscillation maximum. For the normal (inverted) hierarchy, 
the magnitude of the — > v e transition probability is enhanced (suppressed) by about 
10% at Kamioka, 24% at Oki Island, and 37% at L ~ 1000 km in Korea, around the first 
oscillation maximum, |A 13 | ~ tt. When L/E u is fixed at |A 13 | ~ tt, the difference between 
the two hierarchy cases grows with L, because the matter effect grows with E u . Within 
the allowed range of the model parameters, the difference of the A e between SK and a far 
detector at Oki or Korea becomes 



A e peak (L = 653km) - ^ eaj£ (L = 295km) ~ ±0.13. 
A e peak (L~ 1000km) -A; e3k (L = 295km) ~ ±0.26 



(24a) 
(24b) 



where the upper sign corresponds to the normal, and the lower sign for the inverted hi- 
erarchy. The hierarchy pattern can hence be determined by comparing P v . Ue near the 
oscillation maximum |A 13 | ~ tt at two vastly different baseline lengths [26]- [30], indepen- 
dently of the sign and magnitude of sin5 MNS . 

In eq. (I23bl) . it is also found that the first term in B e , which shifts the oscillation phase 
from |A 13 | to |A 13 ± B e \ = |A 13 | ± B e , is also sensitive to the mass hierarchy pattern. As 



in the case for A e , the difference in B e between SK and a far detectors is found to be 
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(25a) 
(25b) 



where the upper sign is for the normal, and the lower sign for the inverted hierarchy. 
This implies that the mass hierarchy pattern can also be discriminated by measuring the 
neutrino energy of the first oscillation maximum at two different baseline lengths. As in 
the case of the oscillation amplitude above, this determination can be made independent 
of the unknown cos5 MNS terms since they cancel in the differences, eq. fl25|) . 

Once the sign of A 13 is fixed by the terms linear in p, the terms linear in A 12 , which 
appear as those proportional to |Ai3| ~ 30Ai2 in eq. (123]) . allow us to constrain sin 5mns via 
the amplitude A e , and cos<5 MNS via the phase shift B e . Therefore, 5 MNS can be measured 
uniquely once the mass hierarchy pattern is determined. From the above discussions, we 
understand qualitatively why the mass hierarchy as well as both sin 5 MNS and cos 5 MNS can 
be determined uniquely by observing the — > v e oscillation probability around the first 
oscillation maximum at two vastly different baseline lengths. Therefore, in order to take 
advantage of this very efficient mechanism to determine all the main unknowns of the three 
neutrino model, one should arrange for high neutrino flux both at near and far detectors 
around the first oscillation maximum which appears at the same energy-to-baseline-length 
ratio, 
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for |5m^ 3 | = 2.35 x 10 _3 eV 2 ; eq. (Ilbj) . This observation led to the T2KK proposal of 
choosing the 3° off-axis beam (£P eak ~ 0.55GeV) at SK and ~ 0.5° off-axis beam (£P cak ~ 
l.lGeV) at L ~ 1000km [261 I2Z|- Since the baseline length of L ~ 653km to Oki Island is 
about a factor of two longer than L ~ 295km for T2K, we may also expect that a beam 
with smaller off-axis angle at Oki Island will enhance its physics capability; see, however, 
discussions in section [61 

The oscillation probabilities for the anti- neutrino, P(u a — > Vp), are obtained from the 
expressions for P(v a — > vp) by reversing the sign of the matter effect term, (p — > —p), and 
that of 5 MNS , (5 MNS y — S MNS ). The differences in the shift terms A e and B e for z/^ — > u e 
and A e and B e for — > v e oscillation probabilities are, respectively, 

A e -A e = 
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where the upper sign should be taken for the normal and the lower sign for the inverted 
hierarchy. From eq. (|27aj) . the difference of the amplitude between v e and v e appearance 
probability grows with the baseline length with the opposite sign for the normal and in- 
verted hierarchies. The dependence on sin 6 MNS is also strong since A e and A e changes in the 
opposite direction when we vary sin# MNS . The phase-shift term B e and B e also change sign 
in the opposite direction for the normal and inverted hierarchies, whose sign is independent 
of cos 5 MNS because of the smallness of the cos S MNS dependence in eq. (127bl) . By the same 
token, the use of both and beams does not improve significantly the measurement 
of cos5 MN g. Because these sign changes occur independently of the /^-dependence of the 
oscillation probabilities, the physics potential of the T2KO experiment will be improved 
by dividing the total exposure time to neutrino and anti-neutrino beams just as in T2KK 



It is worth noting here that all the formalism presented in this section should be useful 
for studying physics potential of T2K plus NOvA [M] , whose baseline length of L ~ 810 km 
lies between those of Tokai to Oki Island (L ~ 653 km) and to Korea (L ^ 1000 km). 



In this section, we study the matter profile along the baselines between J-PARC and 
Kamioka, Oki Island, and Korea by referring to recent seismological measurements [35] - 
[12]. We also show the relation between the off-axis angles at SK and Oki Island, and the 
beam profiles for the relevant off-axis angles [43J . 

3.1 Matter profile along the baselines 

Because Oki Island is placed just at the middle point between SK and the east shore of 
Korea along the T2K beam, as shown in Fig. [H the cross section view of the T2K, T2KO, 
and T2KK experiments along the baselines can be shown on one frame as in Fig. [2j The 
horizontal axis of Fig. [2] gives the distance from J-PARC along the arc of the earth surface 
and the vertical axis denotes the depth below the sea level. Three curves show the baselines 
from the J-PARC to SK (L = 295Km), Oki Island (L = 653km) and to the Korean east 
shore at L = 1000km. The numbers in the white boxes are the average matter density of 
the layer in units of g/cm 3 . 

If we assume that the matter density in each layer of the earth crust as shown in Fig. [2] 
has the value equals to the quoted average, we can estimate the average matter density p 
along the three baselines: 
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Figure 2: The cross section view of the T2K, T2KO, and T2KK experiments along the 
baselines, which are shown by the three curves. The horizontal scale gives the distance 
from J-PARC along the arc of the earth surface and the vertical scale measures the depth 
of the baseline below the sea level. The numbers in the white boxes are the average matter 
density in units of g/cm 3 [35]- [42]. 

The error of these average density can be estimated from the uncertainty of the matter 
density in each region and the location of the boundary of each layer. They are measured by 
using the velocity of the seismic wave in most geophysical researches. From the uncertainty 
in the correlation between the matter density and the measured p-wave sound velocity 
[4"2| [30] . we adopt 6% overall error for the matter density. Small fluctuation in the matter 
density in each layer of the earth crust does not affect the neutrino oscillation probabilities 
significantly, because the contribution from the higher Fourier modes of the matter density 
distribution is strongly suppressed [HI [30]. The locations of the boundaries are measured 
rather accurately from the reflection point of the seismic wave. From Fig. (2J we find that 
the neutrino beam for the Tokai-to-Oki baseline goes through the upper crust layer with 
p = 2.8g/cm 3 , except when it crosses Fossa Magna filled with sediment. The uncertainty of 
the boundary depth between the sediment layer and the upper crust is only a few hundred 
meters. The error from the uncertainty of the boundary depth can hence be safely neglected 
for the T2KO experiment. The traveling distance in the mantle and the crust depends on 
the depth of the boundary between the lower crust and the upper mantle for the T2KK 
experiment. The average matter density is sensitive to the boundary location for the 
baseline of L ~ 1000 km, because the baseline almost touches the mantle as can be seen 
from Fig. [2j The resulting uncertainty has been estimated and is found to be significantly 
smaller than the 6% overall uncertainty [30]. We therefore use the average matter density 
along each baseline given in eq. ( 1281) to estimate the oscillation probability and assign its 
overall uncertainty of 6% in this study. 
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Figure 3: The relation among the off-axis angles at SK, Oki Island, and at a far detector 
in Korea (L = 1000km, 1100km, and 1200km). The horizontal axis gives the off-axis angle 
at SK and the vertical axis is that at the far detector locations. The purple line shows the 
off-axis angle at Oki Island. The red, blue, and black lines are the smallest off-axis angle 
at L — 1000km, 1100km and 1200km, respectively. The points corresponding to the purple 
circles and black squares are used in our numerical analysis. 

3.2 Off-axis angles and beam profiles 

Figure [3] shows the relation among the off-axis angles of the neutrino beam from J- PARC 
observable at SK, Oki, and a far detector in Korea (L = 1000km, 1100km, and 1200km). 
The horizontal axis gives the off-axis angle at SK and the vertical axis gives the correspond- 
ing off-axis angles at far detector locations. The purple line shows the off-axis angle at Oki 
Island, which grows as that at SK is increased, because both SK and Oki Island is located 
in the east (upper) side of the beam center of the T2K neutrino beam; see Fig. [TJ The red, 
blue, and black lines are the smallest off-axis angle observable in Korea at L = 1000km, 
1100km, and 1200km, respectively. In contrast to the Oki Island case, the off-axis angle 
at Korea decreases as that at SK grows, because Korea is in the west (lower) side of the 
beam center, also as shown in Fig. [TJ The points corresponding to the purple circles and 
black squares in the figure are used in our numerical analysis in the following sections. 

In Fig. HI the beam profiles of the three off-axis beams (OAB) with 2.0°, 2.5°, and 3.0° 
off-axis angles at SK are shown at SK, Oki Island (three purple circles in Fig. [3]) and at 
L = 1000km in Korea (three black squares in Fig. [3]), together with the — > u e and 

- 1 v e oscillation probabilities as functions of the neutrino energy, E v in units of GeV. 
The left column panels are for the focusing beam and the right column ones are for the 
Up focusing beam. The upper three rows, (al)-(a3) and (bl)-(b3), display the beam flux 
for the Up and focusing beam, respectively, at L = 1 km in units of 10 10 /cm 2 /10 21 POT 
(protons on target). The first row, (al) and (bl), gives the beam profile observable at each 
far detector for the 2.0° off-axis beam (OAB) at the SK (red solid line), which gives the 



0.6° OAB at Oki Island (green dashed line), and the 1.5° OAB for a far detector in Korea 
at L — 1000 km (blue dotted line). The second row, (a2) and (b2), is for the 2.5° OAB 
at SK, which gives the 0.9° OAB at Oki Island, and the 1.0° OAB at L = 1000 km. The 
case for the 3.0° OAB at SK, which gives the 1.4° OAB at Oki Island and the 0.5° OAB at 
L = 1000 km, is shown in the third row. It is clearly seen that all fluxes at Oki and Korea 
are harder than those at SK. The hardest flux with E peak ^ 1 GeV is found for the 2.0° at 
SK beam at Oki, and for the 3.0° at SK beam in Korea. The 2.5° at SK beam gives almost 
the same flux shape of ~ 1.0° OAB at Oki and at the L = 1000 km location in Korea, as 
shown in Fig. [TJ This may call upon the possibility of locating detectors both in Oki and 
in Korea, where both detectors observe exactly the same OAB around the off-axis angle of 
1.0°. 

The lower two rows in Fig. H]show the — > v e (z/ M — > v e ) oscillation probability in the 
left (right) column, for the normal hierarchy, (PI) and (P3), and for the inverted hierarchy, 
(P2) and (P4). The oscillation probabilities are calculated for the average matter densities 
of eq. f l28|) . and the three- neutrino- model parameters: 

5m\ 2 = 7.5 x 10~ 5 eV 2 , (29a) 

5m 2 13 = ±2.35 x 10~ 3 eV 2 , (29b) 

sin 2 # ATM = 0.5, (29c) 

sin 2 2# SOL = 0.852, (29d) 

with 

sin 2 2# RCT = 0.08, (30a) 

<Ws = 0.0°. (30b) 



The plus (minus) sign in eq. (129bj) is taken for the normal (inverted) hierarchy, whose 
predictions are shown in the fourth (bottom) row of Fig. HI The red solid line in these 
panels is the oscillation probability at the SK, the green dashed line is for Oki Island, and 
the blue dotted line is for a far detector in Korea, at L = 1000 km. 

From the four plots (P1)~(P4) in the bottom two rows in Fig. HI we find that the first 
oscillation peak at SK appears at E v around 0.4 ~ 0.6 GeV for all the four cases, vs 
and for both hierarchies. Since the OAB peaks at around E v = 0.7, 0.6, and 0.5 GeV for 
2.0°, 2.5°, and 3.0°, respectively, as shown by the red solid curves in the first three rows 
of Fig. HI all the OAB at SK between 2.0° and 3.0° can observe — > v e and — > v e 
oscillations with good efficiency. On the other hand, the oscillation probabilities at Oki 
Island are high at E u around 1.0 ~ 1.4 GeV from the green dashed curves in the bottom 
two rows, where the flux is slightly small for the 3.0° OAB at SK, as shown also by green 
dashed curves in the third row. Finally, at L — 1000 km, the — > v e and z/ M — > v e 
oscillation probabilities are large at around E v ~ 2 GeV, and hence the OAB with 3.0° at 
SK is most favorable [26]- [29]. We will confirm the above observations for the whole region 
of the three neutrino model parameter space in the following sections. 




Figure 4: The beam profiles and the oscillation probabilities as functions of the neutrino 
energy. The left column panels are for the focusing beam and the right ones are for the 
Up focusing beam. The upper 6 panels, (al)-(a3) and (bl)-(b3), show the beam fluxes at 
L = 1 km and the bottom 4 panels, (P1)-(P4), show the — > v e (p^ — > v e ) oscillation 
probability with the input parameters of eqs. ( 1281) and ( l29i) for sin 2 2# RCT = 0.08 and 
5 MNS = 0°. The first row, (al) and (bl) gives the beam profile at each far detector for 
the 2.0° off-axis beam (OAB) at SK, the second row, (a2) and (b2) is for the 2.5° OAB 
at SK, and the third row, (a3) and (b3) is for the 3.0° OAB at SK. The fourth row, (PI) 
and (P3), shows the — > v e and — > v e oscillation probability for the normal hierarchy, 
respectively, and the bottom row, (P2) and (P4), is for the inverted hierarchy. 

4 Typical event numbers at each detector 



In this section, we explain how we estimate event numbers at each far detector. 



4.1 Event calculation 



In order to compare the physics capability of the T2KO experiment and that of T2KK 
experiment with a far detector at various baseline lengths and OAB's, we adopt the same 
conditions as those in Refs. [261 123 [30]: All the detectors at Kamioka, Oki Island, and 
Korea are assumed to have excellent detection and kinematical reconstruction capabilities 
for and v e Charged Current Quasi-Elastic (CCQE) events within the fiducial volumes 
of the 22.5 kton at Kamioka (SK) and 100 kton at Oki Island or in Korea. We use the 
neutrino fluxes of J-PARC beam at various off-axis angles [13] and the CCQE cross sections 
off water target [45J to compute event numbers as functions of (reconstructed) neutrino 
energy with the bin width of SE U = 200 MeV at E v > 400 MeV. The energy bin width 
of 200 MeV is chosen to take account of kinematical reconstruction errors due to Fermi 
motion, resonance production and detector resolutions [29] . We take account of background 
from secondary beams, such as u e , v e and z/ M fluxes for the focusing beam, but we do 
not consider other backgrounds including the single tt° production from the neutral current 
(NC) events [29]. In other words, the results of our studies show what a perfect detector of 
a given fiducial volume can achieve with neutrino beams from J-PARC when it is placed 
along the T2K beam line. Reconstruction efficiency and errors as well as the background 
rejection capabilities depend on specific detector designs and their studies are beyond the 
scope of this paper. 

The number of v a CCQE events from the vp flux in the z/^-focusing beam in the n-th 
energy bin, E™ = E^ + (n — 1) x 8E V < E v < E™ + 5E V , at each site are calculated as 

rEJ)+5E v 

N£(isp -> u a ) = M D N A / ^ V0 {E)P v ^ Va {E)a^{E) dE , (31) 

Jed 

where the suffix "D" denotes the place of the far detector (D=SK, Oki, Kr), v^p is neutrino 
or anti-neutrino flavor {y a ,p = v ni v ei^yn^e): is the fiducial mass of each detector 
(M SK = 22.5 kton and M 0ki = M Kr = 100 kton), N A = 6.017 x 10 23 is the Avogadro 
number, $^ is the vp flux from J-PARC [13] which scales as 1/L 2 , P Vfi ^ tVa is the neutrino 
or anti-neutrino oscillation probability calculated by eq. ( II 5 j) with the average matter 
density of eq. ( 1281) . and <r^ CQE is the CCQE cross sections per nucleon in water for each 
flavor type [15]. We also define the event number for focusing beam as 

rE2+SE v _ 

jr D (vp -> u a ) = M D N A / ^ p (E)P^ Ua (E)a2 CQE (E) dE , (32) 

where $„„ gives the vp flux in the focusing beam |43j. 

Because we only consider the contribution from the secondary neutrino fluxes of the 
primary beam as the background in this analysis, each e- and /i-like event numbers in the 
n-th bin are calculated as 

Kb = E + , (33) 
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Figure 5: Typical event number of the e-like CCQE events for E v (E 9 ) > 0.4 GeV at 
SK (22.5kton) with 3.0° off-axis angle and 2.5 x 10 21 POT exposure: (al) and (a2) are for 
z/ M focusing beam, and (bl) and (b2) are for focusing beam; whereas (al) and (bl) are 
for the normal hierarchy, and (a2) and (b2) are for the inverted hierarchy. The results are 
for sin 2# RCT = 0.08, <5 M ns = 

0°, and Pt2k — 2.60 g/cm 3 and the other input parameters 
of eq. f[2Tlj) . The red solid line in each panel denotes the total event numbers, the purple 
dotted line shows the v e — > v e contribution in (al) and (a2), v e — > v e in (bl) and (b2), the 
blue short dashed line gives the sum of the contributions from — > v e and v e — > v e in (al) 
and (a2), that of — > v e and u e — > v e in (bl) and (b2), and the green dashed line shows 
the total background contribution from the secondary beams. 

for a = e,/i at each far detector "D" (D=SK, Oki, Kr). We also define the event numbers 
with focusing beam as 

K,d = Wfa ~* v a ) + Nl{u p -> u a )} , (34) 

/3=e,n 

for a — e, fi at each far detector. 
4.2 SK 

The typical event numbers of the e-like CCQE events for E v (E p ) > 0.4 GeV at SK 
(22.5 kton) with the 3.0° OAB and 2.5 x 10 21 POT exposure is shown in Fig. EJ (al) and 
(a2) are for the focusing, and (bl) and (b2) are for the focusing beam; (al) and 



(bl) are for the normal, and (a2) and (b2) are for the inverted hierarchy. These results 
are for sin 2 2# RCT = 0.08, <5 MNS = 0° and p T2K = 2.60g/cm 3 the other input parameters of 
eq. ( J29l) . In each panel, the red solid line denotes the total event number, which is the sum 
of the signal events, — > v e for (al) and (a2), or — > v e for (bl) and (b2), and the total 
background from the secondary beams shown by the green dashed line. The purple dotted 
line shows the v e — > v e contribution in (al) and (a2), v e — > v e in (bl) and (b2). The blue 
short dashed line gives the sum of the contributions from — > v e and v e — > v e in (al) and 
(a2), and those from — > v e and v e — > v e in (bl) and (b2). 

In Fig. El both the signal (red solid minus green dash) and the total (red solid) number 
of events peak in the first bin (400 MeV < E v < 600 MeV), because both the 3.0° OAB 
fluxes and the oscillation probabilities are largest in the region; see Fig. HJ (a3), (PI), and 
(P2) for — > z/ e , (b3), (P3), and (P4) for z/^ — > v e . The background levels are higher for 
the — > v e oscillation experiments than the — > v e case for both hierarchies. This is 
partly because of the higher level of the secondary beam fluxes in the focusing beam 
[13] . and partly because the CCQE cross sections are larger for vi than the vt {t = fi,e) 

Since the ratio of the — > v e to the z/ M — > v e event numbers is significantly larger 
for the inverted hierarchy than that for the normal hierarchy case in Fig. El one may be 
tempted to conclude that the neutrino mass hierarchy can be determined by using both 
and focusing beams at T2K. This is not the case since the same trend can be expected 
for sin<5 MNS ~ —1, as can be seen clearly from eq. (I27al) . 



4.3 Oki Island 

Figure |6] shows the typical event numbers for the e-like CCQE event at Oki Island (L = 
653 km), where we suppose to place a 100 kton fiducial volume detector. The results are 
for 2.5 x 10 21 POT at J-PARC for both and z/ M focusing beam at the off-axis angle 
of 1.4°, which corresponds to the 3.0° OAB at SK, as shown in Fig. [31 We use eqs. fl29|) 
and (1301 for physics parameters and eq. f!28bl) for the average matter density to generate 



these event numbers. The left two panels in Fig. El (al) and (a2), are for the focusing 
beam, while the right two panels (bl) and (b2) are for the focusing beam. The top two 
panels, (al) and (bl), are for the normal mass hierarchy, and the bottom two panels, (a2) 
and (b2), are for the inverted mass hierarchy. The line types are the same as in Fig. El 
the red solid lines give the total event numbers, the green dashed lines are the sum of all 
the background events from secondary beams, the blue short dashed lines give the sum of 

— > v e and v e — > v e for the focusing beam, that of z/ M — > v e and v e — > v e for the 
focusing beam. The purple dotted lines shows the contribution from v e — > v e {p e — > v e ) 
for {v^ focusing beam. As in the case of SK shown in Fig. El the background from the 
secondary beam contributions are higher for the focusing beam than those for the 
focusing beam. 

In Fig. El we find that the first oscillation peak appears at around 1.0 GeV at Oki Island 
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Figure 6: The same as Fig. El but for Tokai-to-Oki Island (L = 653km), where the fiducial 
volume of 100 kton and the off-axis angle of 1.4° is assumed for the far detector. The results 
are for sin 26 RCT = 0.08 and 5 MNS = 0°, and 

PT20W — 2.75 g/cm 3 , see eq. (128bj) , 



for all the four cases. The peaks at E v ~ 1 GeV in the event numbers are obtained by the 
convolution of the oscillation probability, whose first peak is located around E v ~ 1.2 GeV 
for the normal or slightly above 1 GeV for the inverted hierarchy, as shown in Fig. H] (Pl)- 
(P4), and the 1.4° OAB fluxes that have a peak at around 0.8 GeV in Fig. H] (a3) and 
(b3), as shown by green dashed lines. The difference between the normal and inverted 
hierarchy is larger than that of the T2K experiment. One can observe the second peak 
in the 400^600 MeV bin for the normal hierarchy case, (al) and (bl), but not for the 
inverted case. This is because the matter effect to the oscillation phase, the term B e in 
eqs. (122]) and f l23b|) . grows with the baseline length L, and shifts the peaks of the oscillation 



maximum at |A 13 + B e \ ~ tt, 3tt, ■ ■ ■ in the opposite directions; toward higher (lower) 
energies for the normal (inverted) hierarchy. More accurately speaking, the above phase 
shift pattern applies for cos5 M ns ~ 1, as can be read off from eq. f l23bl) . where the matter 
effect diminishes (enhances) the shift for the normal (inverted) hierarchy. The pattern 
reverses for cos5 MNS ~ —1. Likewise, the matter effects on the oscillation amplitudes are 
also clearly seen: we expect more (less) — > v e events than — > v e events around the 
first oscillation peak when the hierarchy is normal (inverted) as can be seen from eq. (I23ap . 
This pattern is enhanced when sin<5 M Ns ~ — 1 whereas it is diminished when sin<5 M Ns ~ 1- 




Figure 7: The same as Fig. [51 but for Korea (L = 1000km), where the fiducial volume of 
100 kton and the off-axis angle of 0.5° is assumed for the far detector. The results are for 
sin 2 26 l RCT = 0.08 and S MNS = 0°, and p T 20ki = 2.9 g/cm 3 , see eq. (I28cl) . 

4.4 Korea 

In Figs, d we show the typical event numbers for e-like CCQE event at a far detector in 
Korea (L = 1000 km) with a 100 kton fiducial volume detector. The off-axis angle there 
is chosen to be 0.5°, which optimized the T2KK performance with the 3.0° OAB at SK 
[261 1271 129] . The results are for sin 2 29 RCT = 0.8 and S MNS = 0° and the parameters of 
eqs. ( 1291 and ( 1301) . with the average matter density of 2.9 g/cm 3 eq. (128cjl . The panels and 
the lines types are the same as those in Figs. |5] and [6j 

The first peak of the event appears at around 1.6 GeV for the beam for both hierar- 
chies, but that is around 1.8 (1.5) GeV for the z/ M beam for the normal (inverted) hierarchy. 
This is because of the matter effect contribution to the oscillation phase, B e in eq. ( I23bl ). 
and B e obtained from eq. (I23bl) by reversing the sign of p. The magnitude of B e is small 
for cos<5 MNS = 1 in eq. (123bj) . whereas that of B e is enhanced. As expected, the sign of the 



difference B e — B e in eq. (I27bl) reflects the neutrino mass hierarchy. 



Likewise, the difference in the heights of the first oscillation peak is more distinct in 
Fig. [7] for L = 1000 km then that in Fig. [6] for L = 653 km. As expected, the oscillation 
amplitude for the — > v e transition is bigger (smaller) than that for the — > v e transition 
when the hierarchy is normal (inverted). Although the magnitudes of those enhancement 
or suppression factor depends on sin^MNS; can be seen from eqs. ( I23ap and ( I27ap . the 



difference between different baseline lengths among Figs. [51 [61 and [7] depend solely on the 
mass hierarchy pattern, as expressed in eq. fT2"4"l) . 



5 Analysis method 

We introduce a A% 2 function 

= Xstat ± Xsys ± Xpara ) (35) 

in order to compare the physics potential of T2K, T2KO, and T2KK experiments quanti- 
tatively on the same footing. The first term of eq. ( |35|) gives statistical constraints on the 
model parameters from the number of the CCQE events in each bin at each detector: 



Xstat 




,(36) 



where iV" D and iV" D denotes the number of \i- and e-like events, respectively, for the 
focusing beam calculated with eq. (133]) in the n-th bin at each detector, D=SK, Oki, Kr, 
whereas iV™ D and N^ D are for the focusing beam. Their square-roots give statistical 
errors. The summation is over all bins from 0.4 GeV to 5.0 GeV for /i-like events at all 
sites, and for e-like events from 0.4 GeV to 1.2 GeV at SK, from 0.4 GeV to 2.4 GeV at 
Oki, from 0.4 GeV to 2.8 GeV at a far detector in Korea, respectively. The input event 
numbers, (-/V™ D ) input , are generated with eq. fl29|) for both hierarchies, and for various value 
of sin 2 20 RCT and S MNS . We use the average matter density along each baseline of eq. f[2"8"j) 
when we calculate the input event numbers. 

The event numbers in the fit, (-/V™ D ) fit and (iV™ D ) fit , are obtained by varying all the 6 
parameters of the three-neutrino model freely and also by allowing for the systematic errors. 
We consider the following systematic errors in this analysis. We assign 6% uncertainty to 
the overall matter density along each baseline, 

Pd = ffpB , ff = 1-00 ± 0.06 , (for D=SK, Oki, Kr). (37) 

Although we expect positive correlation among f®, we allow them to vary independently 
as a conservative estimate. We assign 3% uncertainty in the flux normalization of and 
focusing beam as 

= /£ = 1.00 ±0.03, (for v a = v^v e ^v e ), (38a) 

(*£.)* = > j£ = 1.00 ±0.03, (for u a = u e , u„ u e ) , (38b) 

respectively, where D=SK, Oki, Kr. We also ignore possible correlations among the flux 
errors. For the CCQE cross sections of neutrino and anti-neutrinos, we assume common 



3% error as 

(a^E))** = f t (a^ E (£)) , ft = 1.00 ± 0.03 , (39a) 

(a£ QE (£)) fit = h (^ E (E)) , h = 1-00 ± 0.03 , (39b) 

for neutrino and anti-neutrino events, independently, but we take f u = f Ve = ft and 
fuu = fv e = ft because of the e-\i universality. Here also, we neglect the correlation 
between Uf and i>t cross section errors. The systematic error for the fiducial volume of each 
far detector is also assigned as 

= f D M D , f D = 1.00 ± 0.03 (for D = SK, Oki, Kr) . (40) 

Summing up, we take account of 32 systematic uncertainties in terms of which Xsys i n 
eq. fl35l) is expressed as 



x 2 



sys 




(41) 



There are 32 normalization factors since for each detector "D" the fiducial volume (/d), 
the average matter density (/?), and the 4 fluxes each for (f® a ) and (f^ a ) focusing 
beam are accounted for, in addition to the overall theoretical uncertainties for vg (fe) and 
vt (ft) CCQE cross sections. 

The last term of eq. fl35|) . Xpara; accounts for the external constraints on the model 
parameters: 

2 / 7.5 x 10- 5 eV 2 - (£mf 2 ) fit \ 2 / Q.852 - sin 2 2flf£ L 

Xpara ~ V 0.2 x lO" 5 J + V 0!025 

| ^ sin 2 2Cc P T Ut -sin 2 2e T y ^ (42) 

The first two terms are from the present constraints from the KamLAND experiment [9] 
listed in eq. ([2]). In the last term, we assume that the new reactor experiments [Ml [T7J 123] 
will measure sin 2 28 RCT with an uncertainty of 0.01 in the near future. We do not impose 
the present constraints on |5^ 2 3 | and sin 2 29 ATM given in eq. ([I]), since the experiments 
studied in this report will measure them more accurately. 



6 mass hierarchy 



In this section, we show physics capability of the T2KO experiment to determine the mass 
hierarchy and compare it with that of the T2KK [2S]-[Sn] and the T2K 12 2 experiment. 
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Figure 8: Minimum Ax 2 of the T2KK , T2KO, and T2K 12 2 experiment to exclude the 
wrong mass hierarchy when only the focusing beam of 5.0 x 10 21 POT is used. Four 
columns give results for <5 MNS = 0°, 90°, 180°, and 270°, respectively, and the horizontal 
axis in each column gives the off-axis angle at SK. The solid (open) square, circle, and tri- 
angle denotes, respectively, the T2KK, T2KO and T2K 12 2 results for the normal (inverted) 



hierarchy. The results are for sin 2 29 h 



0.08 and the other input model parameters in 



eq. (1291) . as well as the average matter density of eq. (1281) along the three baselines. 



Here by T2K122, we examine the option where the additional 100 kton detector is placed 
in the Kamioka site to make the total fiducial volume 122 kton at L — 295 km, which may 
be regarded as a small scale version of Hyper-Kamiokande [16] . 

In Fig. [8], we show the minimum A% 2 of T2KK, T2KO, and T2K122 experiment for the 
Vfj, focusing beam with 5.0 x 10 21 POT. Four columns of Fig. [8] give results for 5 MNS = 0°, 
90°, 180°, and 270°, respectively, and the horizontal axis in each column gives the off- 
axis angle at the SK. The solid (open) square, circle, and triangle denotes, respectively, 
the sensitivity of the T2KK, T2KO, and T2K122 experiments for the normal (inverted) 
hierarchy. The results are for sin 2 29 RCT = 0.08 and the other input model parameters in 
eq. ( 1291) . as well as the average matter density of eq. (1251) along the three baselines. 

From Fig. [HJ we can tell that the physics potential for the mass hierarchy determination 
of the T2KK experiment is far better than the other experiments when the combination of 
3.0° OAB at SK and 0.5° OAB at a far detector in Korea is taken, where the mass hierarchy 
can be determined by more than 4<r level for all the 8 cases (4 values of 5 MNS and both 
hierarchies). The sensitivity of the T2KK on the mass hierarchy reduces significantly as 
the off-axis angle at SK is reduced. This is because OAB with small off-axis angle cannot 
reach Korea when off-axis angle at SK is below ~ 2.5° [26j [27] ; see Fig. [3j 
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Figure 9: The same as Fig. |8l but with both and focusing beams each with 2.5 x 10 
POT. 
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For the T2KO experiment, the mass hierarchy can be determined at the 2a level or 
higher. This capability does not depend strongly on the off-axis angle at SK, because the 
beam intensity around the first peak (E v ~ lGeV) at Oki Island does not change much 
with the OAB at SK, see Fig. SJ 

Furthermore, it is clearly shown in Fig. |8] that the T2K 12 2 experiment does not have 
any sensitivity to the neutrino mass hierarchy pattern for any combinations of 5 MNS and 
the mass hierarchy. This is essentially because the small differences in the oscillation 
probabilities between the normal and the inverted hierarchy can easily be compensated by 
small shifts in the model parameters, such as |^m^ 3 |, sin 2 6' ATM , and S MNS . 

In Fig. [91 we show the minimum Ax 2 of the T2KK, T2KO, and T2K 12 2 experiment to 
exclude the wrong mass hierarchy when both z/ M and focusing beams are used, each with 
2.5 x 10 21 POT to keep the total exposure the same. 

As for the T2KK experiment, shown by the red squares, the improvement is significant 
especially when 5 MNS ~ 0° and 90°, making Ax^^ greater than 20 for all the 8 combina- 
tions of S MNS and sgn(£m 2 3 ), not only for the 3.0° OAB but also for the 2.5° OAB at SK 
|30j . Likewise for the T2KO experiment, shown by green circles, the improvement is most 
significant at 6 MNS = 90° where the smallest grows from ~ 4 in Fig. [8] to ~ 12 in 

Fig. |9j This is essentially because the sin S MNS term in the oscillation amplitude shift A e in 
eq. (I23ap changes sign for the — > v e oscillation as shown in eq. f)27ap for the difference 
A e — A e . This helps resolving the entanglement between the sin<5 MNS dependence ad the 
matter effect in the oscillation amplitudes. Significant improvements are also found for 
5 MNS = 270° by the same reason. Moderate improvements are found for the S MNS = 0° case, 
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Figure 10: The Ax^ in contour plot for the T2KO experiment to exclude the wrong mass 
hierarchy in the plane of sin 2 2# HCT and S MNS . The left figure is for the normal hierarchy 
and the right one is for the inverted hierarchy. The OAB combination for both figures is 
3.0° at SK and 1.4° at Oki Island with 2.5 x 10 21 POT for both and focusing beams. 
Contours for Ax 2 ^ = 4, 9, 16, 25, 36, 49 are shown. All the input parameters other than 
sin 2 26 l R cT and 5 MNS are shown in eqs. ff28|) and ff29l . 

allowing the T2KO experiment with half-and-half and beams to resolve the mass 
hierarchy at 3cr level for the worst case (£ M ns — and normal hierarchy). No significant 
dependence on the OAB at SK is found. 

The capability of determining the mass hierarchy pattern by the T2K122 experiment 
does not appear even by using both and beams. The value of Ax 2 ^ stays almost 
zero, except for just two special combinations; 5 M ns = 270° for the normal hierarchy and 
5 MNS = 90° for the inverted hierarchy. These are the two cases, where the difference between 
the z/ M — > v e and —> u e oscillation amplitudes, A e — A e in eq. (T27aj) is largest or smallest, 
respectively, such that variation of the other model parameters cannot account for the 
difference if the wrong mass hierarchy is assumed. With the same token, Ax 2 ^ exceeds 20 
for T2KO or 30 for T2KK (with ^2.5° OAB at SK) for these two particular combinations. 

So far, we have shown results for sin 2 29 RCT = 0.8 and 4 representative values of <5 M ns; 
0°, 180°, and ±90°. In Fig. UM we show the contour plot of the Ax^in for the T2K0 
experiment to exclude the wrong mass hierarchy in the whole plane of sin 2 29 nCT and 5 MNS : 
The left figure (a) is for the normal hierarchy, whereas the right one (b) is for the inverted 
mass hierarchy. Since no strong dependence of Ax 2 ^ on the OAB angle is found for T2KO 
potential in Fig. [101 we choose 3.0° OAB at SK, which gives 1.4° OAB at Oki, as shown 
in Fig. [3j The results are for 2.5 x 10 21 POT each for and focusing beam. Contours 
in each figure are for Ax 2 ^ = 4, 9, 16, 25, 36, and 49. The input parameters other than 
sin 2 26' RCT and <5 MNS are shown in eqs. ( 128|) and ( 1291) . exactly the same as those adopted 
in Fig. [91 Accordingly the Ax^in values at sin 2 29 RCT = 0.08 agree exactly with those 
presented in Fig. [9] for T2KO with 3.0° OAB at SK, for the 4 representative 5 MNS values. 

It is clearly seen from Fig. [10] that the mass hierarchy pattern can be distinguished at 
3a if sin 2 26'RCT * ^0.09 for any value of i^mns* and for both hierarchies. The most difficult 
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Figure 11: The same as Fig. (TOj but for T2KK experiment with the optimum OAB 
combination, 3.0° OAB at SK and 0.5° OAB at L = 1000km. Ax^ in values are given along 
the contours. 



case is found for 5 MNS ~ 0° for the normal hierarchy, while 5 MNS ~ —135° for the inverted 
hierarchy. On the other hand, the discrimination is easiest at <5 MNS ~ —90° for the normal 
and at 5 MNS ~ 90° for the inverted hierarchy, in accordance with the argument presented 
above for Fig. [9J 

In addition, the contour plots Fig. [10] identify another case at <5 MNS ~ 60° (~ —60°) 
for the normal (inverted) hierarchy, where the difference between the right and the wrong 
hierarchy is large, giving high Ax^ in . The spikes of the contours around these <5 MNS values 
appear as a consequence of the conspiracy between the mass hierarchy dependences in the 
oscillation amplitude A e and in the phase shift term B e . When ^mns* — 60°, the — > v e 
oscillation amplitude shift A e cancels between the matter effect term and the sin <5 MNS term 
for the normal (inverted) hierarchy at around the Tokai-to-Oki baseline; see eq. (123al) at 
L = 653 km. The cancellation is not significant at L = 295 km for T2K, and the two effects 
add up constructively for — > v e oscillation. When the wrong hierarchy is assumed the 
best fit is found for sin 5^ ~ 0, for which there are two solutions, cos 5^ ~ 1 and 
cos 5^ NS ~ — 1, which have significantly different phase-shift; see eq. (I23bp . We find that 
the spike around 5 MNS ~ ±60° in Fig. [H)] occurs when the Ax^ in solution of the wrong 
hierarchy model jumps from cos <5^ NS ~ 1 to cos 5^ ns ~ —1. 

Figure [TTJ shows the same contour plots as Fig. [10J but for T2KK with 3.0° OAB at 
SK and 0.5° OAB at L = 1000km. Significant increase in the Ax^ in values at T2KK is 
clearly seen against those in Fig. [10] for T2KO. Now, the wrong mass hierarchy can be 
excluded at 5a level for sin 2 29 RCT > 0.08 (0.09) if the mass hierarchy is normal (inverted). 
Because the measurement error is dominated by statistics, we find that the 3a sensitivity 
of T2KO with 100 kton detector can be archived with a 40 kton detector for T2KK. 
The CP phase dependence of the T2KK contours is much weaker than that of the T2KO 
contours, especially for smaller sin 2 2# RCT . This is simply because the matter effect terms 
at L <; 1000 km dominate over the sin 5mns and cos 5mns terms, in the correction terms A e 
and B e ; see eq. (j23|) . 
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Figure 12: The Ax^ in contour plot for the T2KO experiment in the plane of sin 2 2# RCT and 
5 MNS when the mass hierarchy is assumed to be normal (left) or inverted (right). Allowed 
regions in the plane of sin 2 2# RCT and <5 MNS are shown for the combination of 3.0° OAB at 
SK and 1.4° at Oki Island with 2.5 x 10 21 POT each for and focusing beams. The 
input values of sin 2 20 RCT is 0.04, 0.08, and 0.12 and 5 MNS is 0°, 90°, 180°, and 270°. The 
other input parameters are given in eqs. fl28l) and ff29l . The dotted-lines, dashed-lines, 
and solid-lines show ^Xmin = 1, 4, and 9 respectively. The shaded region has "mirror" 
solutions for the wrong mass hierarchy giving Ax 2 ^ < 9. 

7 CP phase 

In this section, we investigate the measurement of CP phase S MNS in the T2KO experiment, 
as compared to the T2KK and T2K 12 2 options. In all the cases we adopt 3.0° OAB at 
SK, which makes the SK contribution to the measurements exactly the same, and assume 
2.5 x 10 21 POT each for and focusing beam. 

We show in Fig. [12] the ^-Xmm contour plots for the T2KO experiment in the plane of 
sin 2 28 RCT and 5 MNS , when the mass hierarchy is assumed to be normal (left) or inverted 
(right). The 12 cases are shown in each figure for 5 MNS = 0°, 90°, 180°, 270°, and for 
sin 2 28 RCT = 0.04, 0.08, and 0.12. The other input parameters are given in eqs. (T2"8"|) and 
(|29|) . The allowed regions in the plane of sin 2 26* RCT and 5 MNS are inside of the dotted-, 
dashed-, and solid-contours at A^ 2 ^ = 1, 4, and 9, respectively. The shaded region has 
"mirror" solutions for the wrong mass hierarchy giving A^ 2 ^ < 9, as shown by the red 
solid contours in Fig. [101 Since the sin 2 29 RCT = 0.04 input cases are no longer relevant 
after the measurements eqs. ([7]) and ([8]) by DayaBay [17] and RENO [23], respectively, the 
only parameter regions where we should worry about the mirror solution with the wrong 
hierarchy are around 5 MNS ~ 0° at sin 2 26* RCT ~ 0.08 for the normal hierarchy, and around 
S MNS ~ —135° at sin 2 29 RCT ~ 0.08 for the inverted hierarchy. Since these regions are near 
the 3a boundary, the mirror solutions may be excluded by extending the experimental pe- 
riod or by enhancing the beam power. The sin 2 2# RCT = 0.04 input cases are kept, since they 
show the independence of the <5 MNS measurement error on sin 2 28 nCT clearly. As explained 
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Figure 13: The same as Fig. HH but for T2KK experiment with 3.0° OAB at SK and 0.5° 
OAB at L = 1000km. 



in Ref . [26l [27] . this independence is a consequence of the l/vsin 2 20 RCT enhancement of 
the sin5 MNS and cos5 MNS dependencies in A e and B e , respectively, in eqs. (123aj) and ( 12 3b p . 



which cancels precisely the statistical error which is proportional to a/ sin 2 26 RCT , or the 
square-root of the — > v e and — > u e event numbers. 

It is clearly seen that 5 MNS can be measured with ±20° error for all the 24 cases presented 
in Fig. 1121(a) and (b). This is essentially because the magnitude of the coefficient of sin5 MNS 
in the amplitude shift in eq. (I23a|) and that of cos 5 MNS in the phase shift in eq. ( I23bj) are 



approximately equal. It should be noted that the uncertainty in the sin 2 29 RCT is dictated 
by the external constraint with the error of ±0.01 on the \ 2 function eq. (|35|) . Because 
of the nearly zero correlation between the errors of S MNS and sin 2 29 nCT in Fig. [121 further 
improvements in the precise measurements of sin 2 29 RCT will not reduce the errors of S MNS 
significantly. 

In order to compare the sensitivity of the 5 MNS measurement between the T2KO and 
the T2KK experiments, we show in Fig. [13] the same contours for the T2KK experiment 
with 3.0° OAB at SK and 0.5° OAB at L = 1000km. 

It is clearly seen from the 12 sets of contours, each for normal (left) and inverted (right) 
hierarchy, that the expected error of <5 M ns is ±20° for all the combinations, just as for the 
T2KO experiment shown in Fig. [12] This is remarkable since the event number at a far 
detector in Korea at L ~ 1000 km is significantly smaller than that in Oki at L = 653 km 
because of the flux which decreases as 1/L 2 at long distances. This decrease the overall 
flux is compensated by the wide-band structure of the 0.5° OAB as shown by blue dotted 
lines in Fig. H] (a3) and (b3), which enables the far detector in Korea to observe not only 
the first oscillation peak but also the second one as in Fig. [7J Around the second peak, 
|A 13 /7r| ~ 3, and the sensitivity to sin5 MNS and cos<5 MNS can be three times higher than 
the first peak with |Ai 3 /7r| ~ 1; see eqs. ( 1231) and ( T27j) . In addition, the extended energy 
range covered by the 0.5° OAB allows the T2KK experiment to measure the phase-shifts 
B e and B e accurately, and hence cos£ MNS ; see eqs. ( 123bj) and (127bj) . Indeed, we notice in 
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Figure 14: The Ax^ in contour plot for the T2K i2 2 experiment in the plane of sin 2 2# RCT 
and S MNS when the mass hierarchy is assumed to be normal (ml — ml > 0). Allowed regions 
in the plane of sin 2 28 RCT and <5 MNS are shown for experiments with 2.5 x 10 21 POT each for 
and focusing beam at 3.0° off-axis angle. The input values of sin 2 2(9 RCT are 0.04, 0.08, 
and 0.12 and 5 M ns are 0° (a), 90° (b), 180° (c), and 270° (d). The other input parameters 
are listed in eqs. f[29]) and f l28|) . The red dotted-lines, dashed-lines, and solid-lines show 
AXmin = 1) 4, and 9 contours, respectively, when the right mass hierarchy is assumed in the 
fit, whereas the blue contours give A^ 2 ^ measured from the local minimum value (shown 
besides the x symbol) at the cross point when the wrong hierarchy is assumed in the fit. 



Fig. [TS] (a) and (b) that the error of <5 MNS can be as small as ±15° when 5 MNS = 0° or 180°. 

Finally in Figs. [T4] and US] we show the allowed regions by the T2K 12 2 experiment in 
the plane of sin 2 26' RCT and S MNS for the normal (Fig. [H]) and the inverted (Fig. [15]) mass 
hierarchy, when both and focusing beam at 3.0° off-axis angle are used each with 
2.5 x 10 21 POT. The four 5 MNS cases of 0°(a), 90°(b), 180°(c), 270°(d) are examined for the 
three sin 2 29^ = 0.04, 0.08, 0.12, just as in Fig. [TJfor T2KO and Fig. [TJfor T2KK. The 
other input parameters are also the same, taken from eqs. ( 128]) and ( 129]) . The red dotted- 
lines, dashed-lines, and solid-lines show A;^ 2 ^ = 1, 4, 9 contours, respectively, when the 
right mass hierarchy is assumed in the fit, whereas the blue contours gives A^ 2 ^ measured 
from the local minimum at the blue cross point when the wrong hierarchy is assumed in 
the fit. 

The local minimal value of A^ 2 ^ at the blue cross point is given besides the cross 
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Figure 15: The same as Fig. [H] but for the inverted mass hierarchy (m 2 — m\ < 0). 



mark in the corresponding input sin 2 29 RCT column. As has been explained in section [6] 
and shown in Fig. M, the local Ax^m values are significant only for (5mns — 

-90° when the 

hierarchy is normal (Fig. [H](d)), and for 5 MNS ~ 90° when the hierarchy is inverted (Fig. [15] 
(b)). In order to show the location of the input parameters clearly, we show the global 
minimal point by the red solid blob for each input parameter case. The global minimum 
gives Ax^jjj = in our analysis which ignores fluctuations in the mean number of events 
in each bin. 

We first note that <5 MNS can be constrained uniquely around the above two specific 
cases (<5 M ns — —90° for m| — m\ > 0, or <5 MNS ~ 90° for m| — m\ < 0), since not only 
the wrong mass hierarchy assumption gives non-negligible local A^ 2 ^ as given in Fig. [14] 
(d) and Fig. [15] (b), but also the wrong hierarchy assumption favors the right <5 MNS , with 
slightly larger (smaller) fitted sin 2 29 RCT that compensate for the matter effect for the 
normal (inverted) hierarchy. The la error shown by the red dotted contour is rather large, 
however, about ±35°. In all the other cases, the presence of the wrong hierarchy solutions as 
shown by blue contours significantly reduce the capability of measuring <5mns with T2Ki22- 
More importantly, even if we can remove the wrong hierarchy by other experiments such 
as Nova |34j, the next generation reactor neutrino oscillation experiments [4T] . or by an 
atmospheric neutrino observation with a huge detector [4~8l |4~9] . the T2K122 experiment 
with only one baseline length cannot measure 5 MNS with high accuracy when S MNS ~ ±90°, 
or suffers from the 5$L<, = 180° — 5!™=* solution when <5 MNS ~ 0° or 180° as can be seen 



from the separate red contour islands on the (a) and (c) plots in Figs. [14] and [T5l 

Throughout the analysis of this and the previous sections we fix the OAB at SK at 3.0° in 
order to make the SK contributions to all the three experiments T2KO, T2KK and T2K122 
identical. We find that the performance of T2K122 slightly improves if the 2.5° ~ 2.0° OAB 
is adopted instead, mainly because these fluxes are slightly wider (harder) than the 3.0° 
OAB as shown by red solid curves in Figs. |4] (al)-(a3), (bl)-(b3). 

8 summary 

In this paper, we examine physics potential of a one-beam two-detectors neutrino oscillation 
experiment with an additional 100 kton water Cerenkov detector in Oki Island, which is 
located along the T2K beam line at the baseline length of L = 653 km. Together with 
Super-Kamiokande (SK) at L — 295 km, we can measure neutrino oscillations at two 
different energies for the same oscillation phase proportional to L/E. We may call this 
proposal as T2KO (Tokai-to-Kamioka-and-Oki), whose capability has been compared with 
T2KK (Tokai-to-Kamioka-and-Korea) with L ~ 1000 km for the far detector in Korea, and 
also with T2K 12 2 where the same 100 kton detector is placed at the SK site (L = 295 km). 

As shown in Fig. [TJ since the Oki Island is located in the east side of the T2K beam 
center, just like the SK, the off-axis angle at Oki Island increases as that at SK increases, 
as shown in Fig. [3j The off-axis beam (OAB) with 1.4°, 0, 9°, and 0.6° from J-PARC can 
be observed at Oki Island, when the 3.0°, 2.5°, and 2.0° OAB reaches at SK, respectively. 
The neutrino energy of the first oscillation maximum for the Piv^ — > u e ) and Pip^ — > v e ) 
is between 1.0 GeV and 1.5 GeV, which depend on the CP phase and the mass hierarchy 
pattern, at L = 653 km. Since the {py) beam between 0.6° and 1.4° off-axis angles has 
significant intensity around these energies, as shown by green dashed lines in the upper 
six panels in Fig. [4] we expect that the T2KO experiment can be sensitive to the neutrino 
mass hierarchy and the CP phase, just like the T2KK experiment [24], [26j-[30j. 

For a detector of 100 kton fiducial volume and 2.5 x 10 21 POT exposure each for both 
Up and beams, we find that the T2KO experiment can determine the mass hierarchy 
pattern at 3a level if sin 2 26* RCT = 4|[/ e3 | 2 (l — |t/ e 3| 2 ) is larger than 0.09, by observing 
the e-like CCQE (Charged- Current Quasi Elastic) events; see Fig. [TOJ This result does 
not strongly depend on the off-axis angle of the (z/^) beam at Oki Island, because the 
neutrino intensity at the first oscillation maximum does not strongly depend on the off- 
axis angle at Oki Island. The T2KO sensitivity to the mass hierarchy is about 1/3 (in 
^Xmin) °f tlie T2KK experiment with the optimum OAB combination of 3.0° at SK and 
0.5° at a far detector in Korea with the baseline length around 1000 km. This is because 
the factor of two higher sensitivity of the T2KK over the T2KO experiment as shown by 
eqs. ( 1241) and ( 1251) . which should give a factor of 4 in A^ 2 ^ is partially compensated by 
the smaller average flux by a factor of (635 km/1000 km) 2 ~ 0.4 at a far detector in Korea. 
The sensitivity of the mass hierarchy pattern of the T2K122 experiment, where a 100 Kton 



detector is added at the SK location, is almost zero, except around S MNS ~ —90° (90°) for 
the normal (inverted) hierarchy; see Fig. |9j 

The CP phase in the MNS (Maki-Nakagawa-Sakata) lepton flavor mixing matrix [19J, 
5 MNS can be measured with ±20° error for all the four cases at 5 MNS = 0°, ±90°, and 180°, 
almost independent of the sin 2 29 KCT values [26] 127] as long as the neutrino mass hierar- 
chy is determined. This is because sin(5 MNS can be constrained by the difference between 
the magnitudes of the oscillation probabilities P{v^ — > v e ) and P{v^ — > v e ) around the 
oscillation maximum, whereas cos 5 MNS can be determined by the oscillation phase around 
the first oscillation maximum, or the location of the oscillation peak(s); see eqs. (I23a[) and 
(I27aj) . The sensitivity to the CP phase, S MNS , of the T2KO experiment is similar to that of 
the T2KK experiment, mainly because the smallness of the flux at T2KK is compensated 
by its capability to measure the second oscillation peak at |Ai3/7r| ~ 3 when the sensitivity 
to both cos $mns and sin $mns is a factor of 3 higher than that around the first peak; see 
eqs. 

The T2K 122 option, which may be regarded as a first step toward the Hyper-Kamiokande 
[46], cannot generally determine <5 MNS , mainly because it cannot resolve mass hierarchy by 
itself. Only when <5 MNS ~ —90° for the normal hierarchy (see, Fig. [14] (d)) and when 
5 MNS ~ 90° for the inverted hierarchy (see, Fig. [15] (b)), the constraints for both hierarchy 
assumptions overlap, and the CP phase can be determined uniquely. Even if the mass 
hierarchy is determined by other experiments [M], [2]- [19] the sensitivity to S MNS is rather 
poor at T2K122 as shown in Figs. HH and [15] This is essentially because of the parameter 
degeneracy unavoidable in experiments with only one baseline length, such as those between 
<5 MNS and 180° - <5 MNS when sin<5 MNS ~ 0. 
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